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Abstract 

It is a well-known problem that in supersymmetric models there 
are new CP-violating phases which, if unsuppressed, would give a neu- 
tron electric dipole moment 10 2 to 10 3 times the present experimental 
limit. Here we propose that these new phases are suppressed by CP 
invariance, which is broken spontaneously at a high scale and that this 
breaking shows up at low energies only through a universal phase of 
the gaugino masses. It is shown that this can well fit both e and e' of 
the neutral Kaon system. The electric dipole moments of the neutron 
and the electron should be not much below present limits. A model 
incorporating these ideas in a very economical way is presented. 
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In this letter we propose a model of CP violation that solves the small 
phase problem of supersymmetry and has testable low-energy consequences. 
The idea is that CP violation arises spontaneously at a high scale and is 
communicated to the low-energy world through a common phase of the 
gaugino-masses. All low-energy CP violation would be a consequence of 
this one non-vanishing phase angle. The e parameter of the Kaon system 
would arise primarily from the phase of the gluino mass through the box 
diagram [1] shown in Fig. I. As we shall see, to fit e the phase of the 
gaugino mass must be > 3 x 1CT 3 , and the gluino should be relatively light 
(M~ < 500 GeV). This leads to a value of |e'/e| ~ (1 to 3) x 10" 3 (modulo 
hadronic matrix element uncertainties) which arises in the model dominantly 
via the gluino penguin graph of Fig. 2. Electric dipole moments (edm) [2] of 
the neutron (d n ) [3,4] and electron (d e ) [4] would be induced by the one-loop 
diagrams shown in Figs. 3 and 4, which turn out to be not far below the 
present experimental limits. 

The problem that is solved by this idea is the tendency of the neutron edm 
arising from Fig. 3 to come out about a factor of 10 2 to 10 3 too large in models 
with low-energy supersymmetry [3,4]. In SUSY models there are new sources 
of CP violation in the A and B parameters, the \x parameter, and the gaugino 
masses. If CP is explicitly broken there is no reason, in general, why these 
phases should be small. If one assumes, as is natural, that these phases are 
of order unity and that the various as-yet-unobserved superparticles (gluino, 
squarks) have masses around 100 GeV, then one finds that d n ~ 10~ 22 e-cm, 
to be compared with present upper limit of d n < 10~ 25 e-cm. 

One solution to this well-known difficulty is to assume that CP is a spon- 
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taneously broken symmetry. Then CP-violating parameters are finite, cal- 
culable, and, if they arise radiatively, naturally small. This general approach 
to the problem, which is not new, raises two issues. The first is that spon- 
taneous CP violation leads to cosmic domain walls. These can be rendered 
harmless if they are "inflated away". This requires that CP be broken at 
scales larger than the reheating temperature, which argues for the scale of 
spontaneous CP violation to be much higher than Mw- 

The second issue is how the CP violation arising spontaneously at large 
scales is "fed down" to the Kaon system. Several "feeding-down" mecha- 
nisms have been proposed in the literature [5,6,7]. Those suggested in Ref. 
5 and 6 were motivated by the desire to solve the ^-problem (the strong 
CP problem) using spontaneous CP violation and were therefore necessar- 
ily somewhat intricate. In any event, it was shown in Ref. 8 that these 
non-axion approaches to the ^-problem are fraught with difficulties in the 
context of supersymmetry. In our model the ^-problem is solved by the 
Peccei-Quinn mechanism [9], in particular by the KSVZ invisible axion [10], 
and thus our feeding-down mechanism can be much more straightforward 
than the proposals in Ref. 5 and 6. 

The essential idea is that CP is spontaneously broken by the vacuum- 
expectation values (VEVs) of certain gauge-singlet scalar fields, which we 
will call Si, (Si) ^> M w . These VEVs give large complex masses to some 
vector-like fermions (needed anyway to realize the KSVZ invisible axion) 
which are non-singlet under the gauge group. These fermions, which we will 
denote Q + Q c , do not mix with the known quarks and leptons (owing to 
their PQ charges). When Q and Q c are integrated out, the masses of the 
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gauginos that couple to them will acquire a CP-violating phase at one-loop 
which is naturally of order 10~ 2 to 10~ 3 . 

If grand unification is assumed and it is also assumed that the mass of Q 
and Q c comes predominantly from these gauge-singlet contributions of {Si), 
then the phases of the masses of the SU(3)c, SU(2)l and U(l)y gauginos 
will be very nearly equal. In this scenario the only significant CP violating 
phase in the low energy theory is this common gaugino mass-phase, and thus 
all low-energy CP-violation phenomenology is controlled by one parameter. 

A model will now be presented which shows how this idea can be imple- 
mented in a particularly economical way in which the same sector does the 
breaking of both the CP invariance and the Peccei-Quinn symmetry. 

Consider a SUSY SU (5) model in which, in addition to the known matter 
fields, there is a {5} + {5} with Peccei-Quinn charge —1/2. Denote these 
Q + Q c . Coupling to these are two SC/(5)-singlet superfields Si and S 2 , both 
with Peccei-Quinn charge of +1. A third singlet, S3, carries PQ charge of 
—2. The superpotential of this sector is given by 

W = Q C Q (ASi + / 2 S 2 ) + (anS? + ai 2 SiS 2 + a 22 S 2 2 ) S 3 (1) 

where by CP invariance all the parameters are real. Consider the case where 
(Q) = (Q c ) = (S 3 ) = 0. Then F Sl = F S2 = Fq = Fq = are automatically 
satisfied. The breaking of CP can arise as a result of the F Ss = equation, 

anS 2 + ai 2 SiS 2 + a 22 S 2 = , (2) 

which is solved for 

(S 2 ) = k (Si) , k = (-012 ± \ja\ 2 - 4 ail a 22 ) /(2a 22 ) . (3) 
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If a\ 2 ~~ 4ana22 < 0, then there will be a non-trivial relative phase between 
(Si) and (S2), which breaks CP invariance. 

Note that Eq. (3) leaves the magnitude of (Si) , (S 2 ) undetermined. 
It is easy to fix them at the desired PQ scale in several ways. For ex- 
ample, an extra term in the superpotential, (SiSi — MpQ^j X , where Si 
and X are singlet fields with PQ charges of —1 and respectively, would 
lead to (Si) = M 2 PQ from the F x = equation. The soft SUSY term 
m 2 (\Si\ 2 + |Si| 2 ) will minimize the potential for (Si) = (Si^ = M PQ . For 
an, CL22 and ai 2 all of the same order, one will have (S2) being also of order 
Mpq. Hence the same fields that break U(1)pq, namely Si and S2, break CP 
spontaneously. The relative phase, [arg (Si) — arg (S 2 )], is the source of all 
CP violation in the model, while the phase [arg (Si) + arg (S 2 )] is essentially 
the invisible axion. 

The "quarks" Q + Q c that implement the KSVZ axion idea are the means 
of feeding CP violation to the observable low-energy world. The feeding down 
occurs through the (s)quark-loop contribution to the gaugino mass shown in 
Fig. 5. There will be a soft-SUSY breaking term for the squarks of the form 

Vsoft = A/i(Q c QSi) + A 2 f 2 (Q c QS 2 ) + H.c. (4) 

Thus the phase appearing in the squark mass insertion in Fig. 5 is 
arg (Aifi (Si)* + A 2 f 2 (S 2 )*) while that appearing in the quark mass inser- 
tion is arg (/1 (Si) + f 2 (S 2 )). If A x and A 2 were equal these phases would 
cancel and the one-loop contribution to the gaugino mass would be real. 
However, A\ ^ A 2 in general. Even if A x = A 2 at the Planck scale (as is 
expected in supergravity models), they run differently if an 7^ a 22 and would 
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be significantly different at the Peccei-Quinn scale, Mpq, which we assume 
to be between 10 10 and 10 12 GeV. One finds for the phase of the gaugino 
mass 

, a G ( A l -A 2 \\ 2/i/ 2 |fc|sinA | 

Here M\/ 2 is the common gaugino mass at Mqut, «g the gauge coupling 
strength at M GUT and A the phase of k in Eq. (3). 
A number of remarks are now in order: 

(i) If one neglected the effects of the running of the parameters between 
Mqut and Mpq, then the one-loop calculation of arg(M~) given in eq. (5) 
would be manifestly SU(5) invariant, and the phases of the gluino, the Wino 
and the Bino would be all the same. Interestingly, and slightly non-trivially, 
this result remains true to one-loop order in the RGE even when the running 
is taken into account, as explained below. 

In the momentum range Mpq < ji < Mqut, since SU(5) symmetry is 
not exact, the first term in eq. (1) will split into two pieces, a color-triplet 
(fi) part and an S , f/(2)-doublet (L) part: 

W = Q C Q(/ 1 5 1 + f 2 S 2 ) + L C L(/^! + f 2 S 2 ) + .... (6) 

Similarly, the soft SUSY breaking terms of eq. (4) will split into 

Vsoft = A 1 f 1 (n c ns 1 ) + A 2 f 2 n c ns 2 + a'j^ls,) + a'j>(l c ls 2 ) + h.c. (7) 

At and above M GUT , one has fi = f[, ji — f 2 , and A 1 = A[, A 2 = A' 2 . 
From the renormalization group equations for the various parameters of the 
model we find that in the momentum range between Mqut and Mpq, 



It follows from the above that (A 1 - A 2 ) = (A[ - A' 2 ) and /i// 2 = f{/f 2 at 
all scales. Combining with the scaling of gaugino masses, namely, (ctj/Mj) = 
(ckg/Mi/2), we arrive at the result that the phase of all the gauginos are 
identical at the PQ scale even after 377(5) symmetry breaking. They will 
then remain to be equal down to the weak scale. 

(ii) The gaugino phase, 0, can easily be > 3 x 1CT 3 which is what is 
typically required (as will be seen below) to generate e in the K meson system 
from the graph of Fig. 1. However, Eq. (5) shows that large enough to fit 
e requires that the gaugino masses not be too large, a point that is important 
for the expected magnitudes of the neutron and electron edms. Taking aa — 
1/28 and noting that the magnitude of the function in the curly bracket of 
Eq. (5) is less than unity, we see that a phase angle > 3 x 1CT 3 requires 
(A 1 — A 2 )/M 1 /2 £ 2.1. This is both an upper limit on the gaugino mass and 
a lower limit on the A parameter. Solving the RGE for the A parameters 
we found that (A± — A 2 ) ^ 0.7 Aq, where Aq is the universal A parameter at 
the Planck scale. For Aq = 500 GeV, we see that M\/ 2 < 170 GeV, which 
after RGE corrections correspond to M~ < 500 GeV at the weak scale. The 
experimental lower limit on M~ > 150 GeV implies that A > 150 GeV, which 
could have important consequences for the electroweak symmetry breaking. 

(iii) In general complex VEVs can induce phases in the A and B param- 
eters of the ordinary sector at tree-level as noted in Ref. 6. This can be 
avoided if H/(0 i )| 0i=(0i> is real. Since (Q c ) = (Q) = (S 3 ) = (X) = 0, this 
condition is trivially satisfied for Eq. (1) since W evaluated at fa = {fa) 
vanishes identically. 

(iv) There is no reason to expect any other phase than the common 
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gaugino phase to be significantly large at low energy. For example, in the 
KSVZ axion model such as this, where the known quarks and leptons and 
the Higgs superfields Hi and H 2 have vanishing Peccei-Quinn charge, the 
Q c — Q — Si sector is quite separate from the sector of ordinary matter (except 
through their coupling to the gauge/gaugino particles). Thus no one-loop 
diagram involving (Si) contributes to ji or Bji or to the Yukawa couplings of 
the known quarks and leptons. 

We now turn to the evaluation of the CP violating parameters e, e'/e 
and the neutron and the electron electric dipole moments in the model. The 
AS = 2 CP violating effective Hamiltonian is obtained from the gluino box 
graph of Fig. 1 (the SU(2) gaugino box graph is suppressed by a factor of 
~ 30 relative to the gluino box graph and thus is negligible). 



Here M sq is the (common) squark mass, a, f3 are the color indices, <f) is the 
phase of the gluino mass [Eq. (5)], x = M~/M^ q and the function f(x) is 
defined as 



with /(l) = 1. The parameter 5l R is defined to be 5lr = m~ ~ jM\ ' . Since 
the mass-splitting among squarks is constrained phenomenologically to be 
small, we have treated the d^ — s R mass insertion (denoted by m~ ~ ) in Fig. 
1 as small perturbation. Note that other gluino graphs which do not involve 
d~LSR mass insertions (e.g., one with d^s^ mass insertion) are real and do not 
contribute to e. This simplification is a consequence of the fact that only 



-A5=2 
cff 




5 2 LR sm2(p xf(x) -s Ra d a L s Rp d{ + -SRadlsRpdl - (L R). (9) 




(10) 
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the gluino mass has a non-vanishing phase. The contribution to ReMi2 (or 
Amjf) from Fig. 1 is obtained by the interchange sin20 <-> cos20 and taking 
a relative plus sign between the (LR) and (RL) terms. 
The e parameter evaluated from Eq. (9) is given by 




where in the vacuum saturation method of evaluating the K — K matrix 
element B would be 1 by definition, rj is the QCD correction factor from 
M sq to the hadronic scale. If a s in Eq. (11) is evaluated at the \x = M sq: 
then i] ~ 1.8[11] for a s ~ 0.12. fx — 165 MeV is the Kaon decay constant. 
The function xf(x) is slowly varying with its value changing from 1 for x — 1 
to 1.1 for x = 0.1. Using m s = 150 MeV, m d = 10 MeV and x = 1, we 
obtain by fitting |e| = 2.3 x 10" 3 , 

5sin20(5L - 5l L ) ~ 3.2 x 10- 9 (^^7) • (12) 

In our spontaneous-CP violation mechanism, the phase angle is naturally 
of order 3 x 10~ 3 , so that the larger of the mass-splittings, 5lr or 5rl must be 
~ 10~ 3 , assuming that they are not accidentally close in value, (i.e., assuming 
that one of them dominates). Demanding that the contribution from the real 
part of Fig. 1 not be larger than the experimental value of Aro/f, we obtain 

Bcos2<t>{5\ R + 5\ L ) < 1.0 x 10- 6 ( 3QQ ^ ) . (13) 

From Eq. (12) and (13), we obtain the constraint > 3 x 10 -3 . 

The dominant contribution to the AS = 1 CP violating effective Hamil- 
tonian arises from the gluino penguin graph of Fig. 2. (The U(l)y gaugino 
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penguin contribution is two orders of magnitude smaller, the SU (2) gauginos 
do not contribute directly. There are graphs involving W + Hi mixing, but 
these are an order of magnitude smaller.) Evaluating Fig. 2 we obtain 

K A / =1 = Q ^^^-(SLR-SR^sm^gix) [s\ a ia,„(l - l5 )dG% + H.c 

(14) 

where 

, , 2 (2 + 3x - 6x 2 + x 3 + 6x\nx) , , 

glx) = (T^r (15) 

with g(l) = 1. We use the bag model calculation of Ref. (12) to evaluate 
the hadronic matrix element in Eq. (15) and obtain 

t = 5.6 x l^B!rfw^y/ig{x){S LR - 5 RL ) (^^) • (16) 

Here rj' is the QCD correction factor, rj' = [a s (fx) / a s (M)} - 92 ~ 3, where 
pL ~ 1 GeV and a s (/j,) ~ 0.4 has been used. B' is a factor introduced to 
parameterize the uncertainty in the matrix element and is defined to be 1 if 
the Bag model matrix elements given in Ref. (12) are exact. Combining Eq. 
(12) with Eq. (16), we obtain the prediction (for x — 1) 

e' _ _ « (B'\ ( M sq \ 1 



- ~ 2.7 x 10"° — 



0^4^- (17) 



e \B J \300 GeV J 5 LR + 5 RL 

If we make the reasonable assumption that either 5l R or S R l dominates the 
squark mass-splitting, we obtain \e'/e\ ~ (1 to 3) x 10 -3 . This is clearly in the 
range suggested by experiments. Note that the sign of e'/e is not predicted 
in our model. 

In supergravity models, if the minimal supersymmetric spectrum extends 
all the way upto the Planck scale, the squark mass-splitting will be too 
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small for Eq. (11) to account for e. However, the MSSM spectrum is not 
expected to hold all the way to M Ph since the GUT threshold will in general 
bring in new effects [13]. A simple example is the realization of the see-saw 
mechanism for neutrino masses. Between the GUT scale and the Planck 
scale, the Dirac and Majorana neutrino matrices, with their elements not 
necessarily small, will contribute to the evolution of the squark mass matrix. 
The running in this short momentum range can result in relatively large 
values of the mass-splitting. A typical diagram which can generate (IlSr 
mixing via the neutrino Dirac mass matrix in SU(5) is shown in Fig. 6, 
which can lead to Slr ~ (10~ 4 to 10~ 3 ). It has also been emphasized [8] 
that the squark mass-degeneracy in supergravity models, in the absence of 
additional symmetries, will naturally be 5m 2 sq /Ml q ~ 0(a). In realistic string 
compactification scenarios, the squark degeneracy is indeed of this order. 

One of the most interesting consequences of our fundamental hypothesis 
that all low-energy CP violation is the result of a common gaugino-phase is 
that both d n and d e are to be expected at a measurable level. Of course there 
are large hadronic uncertainties in d n , but it is generally estimated that, with 
phases of order unity and sparticle masses of order 100 GeV, d n from Fig. 
3 will be about 10 2 to 10 3 times the experimental bound, as noted earlier. 
Since we require our gluino phase to be > 3 x 10~ 3 to fit e, it is natural to 
expect d n to lie not far below the present bound. 

Again, there are too many presently unknown SUSY parameters involved 
to allow a calculation of the electron edm. However, as noted in Ref. 14, if 
all the superparticles have comparable masses and the gaugino phases are all 
comparable, one would expect that d e ~ 10~ 2 d n . Thus d e should lie not far 
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below 10 e-cm. 
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Figure Captions 

• Fig. 1. A box diagram whereby the phase of the gluino mass con- 
tributes to 7m(Mi 2 ) in the neutral Kaon system. 

• Fig. 2. The gluino penguin graph contribution to \e'/e\. 

• Fig. 3. A contribution to the edm of the <i-quark coming from the 
phase of the gluino mass. A similar diagram exists for the w-quark. 
These in turn induce an edm of the neutron of comparable magnitude. 

• Fig. 4. A contribution to the edm of the electron arising from the 
phase of the photino mass. There are several other diagrams involving 
neutralinos and charginos. 
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• Fig. 5. The diagram by which the gaugino masses acquire a phase of 
order 3 x 10~ 3 . 

• Fig. 6. A diagram contributing to the (IlSr squark mass insertion 
proportional to the neutrino Dirac mass matrix. 
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